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SUMMARY 


Transient and steady state combined natural and forced convective flows over two in-line 
finite thickness fins (louvers) in a vertical channel are numerically solved using two methods. The 
first method of solution is based on the "Simple Arbitrary Lagrangian Eulerian" (SALE) 
technique which incorporates mainly two computational phases: (1) a Lagrangian phase in which 
the velocity field is updated by the effects of all forces, and (b) an Eulerian phase that executes 
all advective fluxes of mass, momentum and energy. The second method of solution uses the 
finite element code entitled FIDAP. In the first part of this study, comparison of the results by 
FIDAP, SALE and available experimental work were done and discussed for steady state forced 
convection over louvered fins. Good agreements were deduced between the three sets of results 
especially for the flow over a single fin. In the second part of the study and in the absence of 
experimental literature, the numerical predictions were extended to the transient transports and 
to the opposing flow where pressure drop is reversed. Results are presented and discussed for 
heat transfer and pressure drop in assisting and opposing mixed convection flows. 


INTRODUCTION 

Louver arrays are used to enhance the performance of compact heat exchangers. If the orientation 
of the exchanger is vertical and the flow rates are low, the buoyancy forces would effect the heat 
transfer and pressure drag characteristics of the fins. Mixed convection near rectangular fins 
with finite thickness has been studied by Kurosaki et al [1], Sparrow et al. [2], and Suzuki et al. 
[3]. Reference [1] has provided experimental data for a single fin, two collinear fins, two parallel 
fins and a staggered array of fins. Suzuki et al. [3] presented finite difference solutions for an 
array of very thin fins in assisting (upward) flow and discussed heat transfer characteristics of 
arrays. Transient mixed convection over a single fin was studied by Khalilollahi and Joshi [4] 
where temperature overshoots and enhanced heat transfer rates were observed for higher Grashof 
numbers. The transient and steady state assisting flow over two in-line was numerically 
investigated in Ref. [5] and steady state results were compared with some experimental data 
reported in Ref. [1]. 


^Assistant Professor of Mechanical Engineering 


PftSCfiDiNG PAGE BLANK NOT FILMED 


127 




Figure 1. Array geometry and computational domain (L=.022 m, 
S=.015 m,H=.105 m,b=.01 m,t=.0022 m) 


The present study is intended to (1) enhance the confidence in the solution by the finite 
difference FORTRAN code, SALE (Simple Arbitrary Lagrangian Eulerian) through comparisons 
with the solution by the finite element package, FIDAP, and (2) to investigate the opposing 
convective flow where the pressure difference between the top and bottom sections reverses 
causing pressure field to oppose the buoyancy force. 


NUMERICAL PREDICTIONS 
Finite Difference Scheme 

Figure 1 shows the model geometry, flow domain and computational grid to the left of x-axis. 
This domain is used in both SALE and FIDAP solutions. The flow field is governed by the 
conservation equations in dimensionless form, 
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dUldX+dVfdY=0 


( 1 ) 


dUfdx * UdUldX + VdUldY = -3P/3X +dUjdx+dR x JdY 

dvjdx * udvjdx+vdvidY = -ap/ay^aiy ar+aiy ax +Gre 
ae/ax +c/ae/ax+yae/ar=[a 2 e/ax 2 +a 2 e/ay z ]// , r 


where 


n II = 23 f//ax , n^^dvjdY , n^acz/ay+atyax , 
Gr=gP(r ff -r c )L 3 /v 2 , Pr-\i-C p lk , Nu=qL/AT k , Q=(T-T C )I(T„-T C ), 
X=xfL, Y=ylL, U=uLjv, F=vl/v, t=rv/L 2 , P=pL 2 lpv z 


The boundary and initial conditions are 

t=0, £/=F=0=O, ( initial conditions) 

T>0, 

a^/ax=c/=ae/ax=o (*r x=o, #l) 

£/=F=0, 0=1 (on fin surfaces) 

[/Vo “ [ p ]r.n =AP (pressure drop) 


The axes of symmetry are at X=0 and X=b/L, as shown in Fig. 1. The above equations were 
solved numerically by a Fortran code which incorporates Simple Arbitrary Lagrangian Eulerian 
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finite difference scheme. This technique is described by Amsden et al. [6]. SALE procedure 
includes a Lagrangian explicit method where computational cells move with the flow, and an 
Eulerian phase in which the cells are returned to the original position. This phase estimates the 
effects of advective fluxes of mass, momentum and energy on the flow parameters. 


Finite Element Scheme 

The finite element solution of conservation equations governing the laminar flow with 
boundary/initial conditions in this problem was made possible through the available educational 
version of the FIDAP package [7]. The grid independence was determined by doubling the 
number of elements until less than 2% difference in maximum velocity at midsection of the lower 
fin was observed. The same procedure was incorporated in the finite difference scheme, SALE. 
An acceleration factor of 0.3 , the quasi-Newton solver, pressure penalty formulation, and a 
fixed time increment (.05 sec) were used. For both schemes, similar unequally sized grid was 
assigned with higher cell density near the heated fin surfaces (Fig. 1). The convergence in all 
cases was relatively fast and smooth. 


RESULTS AND DISCUSSION 


Some transient problems in natural and mixed convection flows have been previously solved by 
using SALE procedure with favorable results [8,9]. This study intends to apply this technique 
and finite element analysis to the phenomenon of heat dissipation of in-line finite thickness 
louvers. In the absence of empirical data, the predictions of transients and thermal characteristics 
of fins in adverse pressure fields can be valuable. In addition. Comparison of the two set of 
predictions (by FIDAP and SALE) can evaluate the reliability and accuracy of these methods 
when applied to thermal design problems. 

Figure 2 indicates the heat transfer steady state performance of a single rectangular fin in 
assisting flow where positive pressure drop assists the buoyancy force. Reynolds number range 
is between 30 to 600. This is common for air-cooled compact heat exchangers. The aspect ratio 
t/L has been shown to have minor effect on the overall heat transfer rates [1]. The results are 
shown for Nusselt number vs. Grashof number and for a fin with t/L = .2 in infinite (air) 
medium. The correlation for forced convection flow over a flat plate is (shown in Fig. 2) 


Nu = .644 Pr 1/3 Re 1/2 


(5) 


The agreement between the three sets of predictions is good especially in the range where 
buoyancy is dominant and Re number is low. The onset of disturbed flow at higher flow rates 
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Figure 2. Steady state Nusselt number for single fin 
in assisting flow 



Figure 3. Comparison of steady state Nusselt number for 

assisting flow 
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and adjacent to fin comers may be accounted for the small disagreement at Re >200. The 
experimental values seem to be average of the two numerically predicted sets. 

Figure 3 presents the comparison of Nusselt number for assisting steady state flow over two 
fins (as shown in Fig. 1) which is positioned in an array. The presence of neighboring fins in 
the array creates the symmetry lines, Y =0 and Y =b/L. The higher values of Gr/Re 2 represents 
the dominance of natural convection while lower range of Gr/Re 2 indicates the forced convection 
regime. The experimental data from Ref. [1] is for in-line fins in infinite medium. This explains 
the difference between the numerical and experimental values for the downstream fin (#2). 
However, predictions of FIDAP and SALE are in favorable agreement especially in buoyancy 
dominant regime. For reference, the relation for mixed convection heat transfer for a flat plate 
is shown in Fig. 3 [10]. 

The estimation of pressure drop is an important consideration in design of heat exchangers. 
In the absence of experimental data for fins of present study, numerical predictions by SALE and 
FIDAP are presented and compared -.for pressure coefficient and for pure forced (Gr=0) and 
mixed (Gr=1480) convections as shown in Fig.4. The agreement is very good especially for 
lower Reynolds number (Re<40) where convection starts by buoyancy action. The difference 
is steady but higher (about 10%) for higher Re numbers where forced convection is dominant. 



Figure 4. Comparison of pressure coefficient for 
assisting flow 
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Figure 6. Comparison of pressure coefficient for 

opposing flow 
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Figure 7. Comparison of transient Nusselt number and average 

velocity for assisting flow 



Figure 8. Comparison of transient Nusselt number and average 

velocity for opposing flow 
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™f K te ™ opposing flow " 1S § lven t0 ±e case where the pressure drop is reversed. That 
means the higher pressure at the top section presses the flow downward while opposing the forces 

tt bU °^ a ?, Cy ' CaSe P resents 111 interesting discontinuity trend for the flow and heat transfer 
Upward free convection begins with AP=0 at about Gr/Re 2 =1.3 (Re=34) causing higiS he^ 

f ° r thC downstr ®f m fin (#1 )- Flow st ays upward but looses intensity when the 
magnitude of pressure drop (which is negative) increases. Diminishing flow corresponds to the 

i°h n R !i 10 Sh T n Fig - 5 ‘ InCreasing ** of AP wilUhSge the 

direction of the flow downward not gradually but stepwise where a sudden shift is observed for 

die Nu number of both fins. Thereafter the flow starts from the upper limit ,Gr/Re 2 = 1 3 in 

downward mixed convection region and moves to lower Gr/Re 2 region approaching Dure forced 

^ differ Tn eS betWeen SALE FIDAP are ^3^ 

sfmfl n a r n ^nd UP T T/ 68 ™ 3,1(1 for both fms in * e forced convection dominant region, 
of f ° r preSSUre coefflcient ^ opposing flow (Fig. 6). Lower magnitude 

w . , ’ ( . f ^ | 1S no * sufficient to push the flow downward in the region where Re < 33 

With increasing lAPl The flow vanishes (Re=0) but then suddenly ch5iges diStfon to a 
downward flow with Re of about 33. Increasing the magnitude of the pressure difference 

O^eSTthP 6 C dlfferenC6S between FIDAP SA L F results become more noticeable 
Overall the agreement seems to be favorable. 

A P „ heat K transfer f nd avera § e verticaI component of velocity for Gr=1480 and 

±3484 are shown in Figures 7 and 8, where predictions by SALE and FIDAP are 

tr„ s the r siem Nu vs ' dimensi ° n,ess •*■«. * . and asking n 0 » 

(AP- +3484) The solutions show minima at early times, then reaching steady state values 
quickly especially for the downstream fin.These minima are caused by the onset of convection 
and after the early conduction heat transfer lowers. The agreement seems to be fairlv good Thp 

0f ? by ™ AP f0r fi " #1 is » the V?ry 

good agreement is observed for the average velocity (V=vL/v). ^ 

Figure 8 represents the same trends observed in Fig. 7. The flow is an opposing tvDe since 
iP is negative Initially the flow is negative or downward, but eventually^ reaches steadv 
upward flow since buoyancy effects are more dominant. A early times ( 05 < 7? 15) uSer 
downward pressure force , die flow is downward before the strengthening of buo™'cy effects 
This region presents higher Nu values for fin #2, since it is die upstream L S7the C 
direction changes while approaching steady state with fin #1 becoming the upstream fin The 
torepancy between FIDAP and SALE is moderately significant for L fettaTd at abou 

( 05 <;<^ 5 i S' average Vel0Ci ' y iS fine exce P' at intermediate times 
(.Ub < t < .25) when FIDAP overestimates the predictions by SALE. 


CONCLUSION 


for TT 16 "* T 1 SKady T* heat ,nmsfer characteristics Of mixed convection were analvzed 
for assisting and opposing flows over two in-line vertical isothermal fins sms,, „ ,* ■ “ 

flow was predicted by three means, namely FIDAP and SALE that are numerical schemes^ 
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bv available experimental data in the literature. The trends in pressure coefficient were presented 
for above cases A discontinuous trend in Nu number, flow rate, and pressure coefficient was 
observed for negative AP values (opposing flow). Overall predictions by the two schemes, 
FDDAP and S ALE , compared favorably. Future work is planned to study further this 
discontinuity and to conduct experiments on verification of numerical results especially for 

opposing flows. 
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